Effect of the Tunneling Conductance on the Coulomb Staircase 
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Quantum fluctuations of the charge in the single electron box are investigated. The rounding of 
the Coulomb staircase caused by virtual electron tunneling is determined by perturbation theory up 
to third order in the tunneling conductance and compared with precise Monte Carlo data computed 
with a new algorithm. The remarkable agreement for large conductance indicates that presently 
available experimental data on Coulomb charging effects in metaUic nanostructures can be well 
explained by finite order perturbative results. 
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Coulomb blockade effects in metallic nanostructures 
are well understood in the region of low temperatures and 
for small tunneling conductance Gt ^ Gk = /h [Q. 
Recently, several groups have started to investigate in de- 
tail the breakdown of charging effects. High-temperature 
anomalies due to weak Coulomb blockade have been stud- 
ied [|-|, and progress in fabrication techniques has lead 
to reliable experimental data for systems with a tunneling 
conductance of several Gk f^-^10|] . Most of the theoreti- 
cal work on Coulomb blockade for strong tunneling can 
roughly be divided into two groups. On the one hand, 
a significant body of work restricts the theory to usually 
two charge states pT|-p^ . This requires the introduction 
of an arbitrary cut-off that enters the final results and 
deranges the comparison with experimental findings. On 
the other hand, systematic perturbative calculations in 
powers of the tunneling Hamiltonian Ht can be shown 
10 to be independent of the electronic bandwidth and 
give results in terms of experimentally measurable quan- 
tities. While this latter approach was successful [ p^jT6| 
in explaining some of the recent data on strong tunneling 
in the single electron transistor the range of valid- 
ity of perturbative results is not known a priori. In this 
Letter, we focus on the simplest system displaying charg- 
ing effects, the single electron box. Perturbative results 
up to sixth order in Ht are derived and compared with 
Monte Carlo data. We shall demonstrate that perturba- 
tion theory indeed does remarkably well. 
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FIG. 1. Circuit diagram of the single electron box, consist- 
ing of a tunnel junction in series with a capacitor. 

The circuit diagram of the single electron box in Fig. |l| 
shows a tunnel junction with capacitance Ct and tun- 
neling conductance Gt in series with a capacitance Cg 
biased by a voltage source Uc^. Provided the tunneling 
conductance Gt ^ Gki the island charge q is quantized, 
i.e. q — —ne, where n is the number of excess electrons 
in the box. At zero temperature, n is a staircase func- 
tion of the external voltage J7ex- This idealized behavior 
is modified at finite temperatures, where the occupation 
of higher charge levels leads to a smearing of the stair- 
case jlTt . Furthermore, the finite tunneling conductance 
causes a hybridization of the lead and island electronic 
states so that the island charge is no longer strictly quan- 
tized. Since the range of validity of perturbative results 
is most restricted at zero temperature, we focus attention 
here on this case. 

At zero temperature, the ground state energy £ as a 
function of riex = C'ct^x/e determines the average island 
charge 
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(n) = ncx - TTTT 
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Here Ec = e^/2C is the classical charging energy, where 
the island capacitance C is the sum of the capacitance 
Ct of the tunnel junction and the gate capacitance Cg, 
cf. Fig. 1^. Formally £ may be written as a perturbation 
series in powers of the tunneling conductance, leading to 
a diagrammatic representation of £ The zeroth or- 
der term of the ground state energy for vanishing electron 
tunneling is determined by the minimum of the electro- 
static energy and reads £q — EdriQ ~ tIcx)^ where uq is 
the integer closest to ricx- Hence, as function of the ap- 
plied voltage, the island charge displays the well known 
Coulomb staircase (n) = uq. Because of the periodicity 
and antisymmetry of (n) as a function of riox, we may 
confine ourselves to < Uex < i- 
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Since the tunneUng Hamiltonian Ht transfers an elec- 
tron from the lead to the island electrode or vice versa, it 
has no diagonal components in the basis of unperturbed 
energy eigenstates and only even orders in Ht contribute 
to the perturbation series of £. The second order term 
gives a contribution proportional to the dimensionless 
conductance a — Gt/Gk and is represented by the two 
diagrams shown in Fig. The arc to the right describes 
the formation of a virtual electron-hole pair in the inter- 
mediate state with an excess electron on the island, that 
is n = 1, and a hole in the lead electrode. The second 
diagram describes the corresponding process with an in- 
termediate state of reduced island charge n = — 1. In the 
interval < ricx < i considered, the two diagrams give 
a contribution to the average island charge of the form 

w 
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where g = a/47r^. The contribution of fourth order in 
Ht follows from the diagrams depicted in Fig. ||b. Here 
the first eight diagrams describe processes with two in- 
termediate electron-hole pairs created and annihilated in 
all possible ways. The remaining diagrams have a lower 
order diagram inserted, as indicated by the prolongation 
of the arc across the vertical line. These diagrams cor- 
respond to terms in the Rayleigh-Schrodinger perturba- 
tion series with energy denominators squared. Each of 
the arcs depicts an integral over the energy of a virtual 
electron-hole pair with a spectral density which becomes 
linear in the infinite bandwidth limit. An element of the 
vertical line corresponds to an energy denominator de- 
scribing the energy difference between the intermediate 
virtual state and the ground state. The diagrammatic 
rules are explained in detail in [p^ and it was shown 
there that although each single diagram diverges in the 
infinite bandwidth limit, the sum of all diagrams of a 
given order remains finite. The contribution of second 
order in a to the average island charge reads Il9| 
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Here s.t.(— ricx) stands for the same sum of terms with 
riox replaced by — riex showing explicitly the asymmetry 
of (n) in the applied voltage C/cx- Further, Li2(a::) = 
— /p" dz \n{l — z)/z denotes the dilogarithm function [ pO[ . 

In third order in a, one has to evaluate 160 diagrams, 
some of which are shown in Fig. |^c. There are 80 di- 
agrams without insertions, such as the left diagram in 
Fig. 64 diagrams with one insertion, and 16 diagrams 



with two insertions. Using the diagrammatic rules, we 
have to deal with three-fold energy integrals of rational 
kernels. These integrals can be done analytically lead- 
ing to polylogarithms pO| and powers of logarithms of 
rational arguments. Since the full analytic expression is 
rather involved , we present explicitly only results for 
two limiting cases rtcx — ^ and ncx — > 
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FIG. 2. Diagrams for the ground state energy in a) first or- 
der, b) second order, and c) selected diagrams in third order 
in a. 

For small external voltage, the average island charge 
grows linearly as 



(Q) =e{n) = C*U,: 



(4) 



where C* is an effective capacitance of the box. In the 
absence of Coulomb blockade effects C* = Cg , while for 
strong Coulomb blockade, i.e., in the limit of vanishing 
tunneling conductance, C* = 0. It is thus natural to 
characterize the strength of the Coulomb blockade effect 
by an effective charging energy E* defined by pa 



Ec C dricx 



d(n) 



The perturbation series gives 

J7'* 

-^^l-4g + V- 



Bg' + Oig' 



(5) 
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where A — 5.066... and B — 1.457... are analytically 
known coefficients whose explicit form is too lengthy to 
present here. 



2 



The perturbation series is well behaved in the region 
of interest, except in the vicinity of riox = where log- 
arithmic divergences arise. This unphysical behavior in- 
dicates the failure of perturbation theory due to the de- 
generacy of the ground state in this limit. For riex — > 5 
one finds from the analytical expression for the average 
island charge 

,2 I l,„3 /„ I /;„2 



(n) = ag^ + bg'' - {g + 6g^ 
-{2g^ + 24g^)ln^d-- 
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where S = ^ — ricx, and where the coefficients a, b, and 
c read numerically a = —9.7726..., b = —70.546..., and 
c = 65.462... . The leading order logarithmic terms in Eq. 

are —g In (5 — 2g^ In^ S — Ag^ In^ 6. These terms come 
from the diagrams shown in Fig. |^, where all intermedi- 
ate states are confined to the two charge states n = 0, 
1 which are degenerate at riex = 5- The most divergent 
logarithmic term of order k reads —^{2gln6)'^, and we 
get by resummation 

— g In S 



{n) = 



1 



2g\n6 

This result was previously derived by Matveev 
renormalization group techniques. 
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FIG. 3. Diagrams of order a'''^^ for the ground state energy 
giving leading order logarithms. 

To explore the range of validity of these higher order 
perturbative results we have carried out precise Quantum 
Monte Carlo (QMC) simulations for the single electron 
box. In contrast to earlier attempts ||T^;^,^ , we do not 
work in the phase representation, but simulate configura- 
tions directly in the charge representation, keeping track 
of all intermediate electron-hole pairs. A general numeric 
scheme for evaluating series of integrals directly in the 
continuum, i.e., without invoking artificial discretization 
of the integration variables was explained in |24|j2^ ] . It is 
possible then to develop an exact (without systematic er- 
rors) QMC algorithm summing diagrammatic series . 
Here, the integration variables are imaginary times of 
charge transfer events, and we have employed this "di- 
agrammatic" QMC for summing all graphs for the par- 
tition function llj]. The efficiency of the new algorithm 



has allowed us to simulate the single electron box at ex- 
tremely low temperatures {f3Ec as large as 10*). While in 
the phase representation ricx ^ results in a sign prob- 
lem, in the charge representation all contributions are 
positive definite even at finite external voltage. This has 
enabled us to obtain for the first time QMC data for the 
entire staircase function. 
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FIG. 4. Effective charging energy as a function of the di- 
mensionless conduction a = Gt/Gk- Perturbative results 
are compared with QMC data. 

In Fig. |[ the effective charging energy E* at zero tem- 
perature is shown as a function of the dimensionless con- 
ductance a. Apart from the point for a = 21, the QMC 
data are consistent with previous results by Wang et al. 
p2| , obtained with a different algorithm. For small ex- 
ternal voltages, the analytical result to second order in a 
is correct with errors below 4% for values of the dimen- 
sionless conductance up to a w 8, while the third order 
result extends to a ~ 16. 

The range of validity of perturbation theory shrinks 
with increasing external voltage and is worst near riox — 
i. In Fig. ^, the rounded Coulomb staircase is depicted 
for a = 5 and a = 10. While the analytical result di- 
verges at riox = 5, we find that third order perturbation 
theory in a remains valid with errors below 4% up to 
ricx ~ 0.495 for dimensionless conductance a — 2 (data 
not shown), up to nex ~ 0.45 for a = b, and up to 
rtcx ~ 0.4 for a = 10. Since for tIox — 0.45 the charg- 
ing energies for n — Q and n — 1 differ only by Q.lEc^ 
deviations from the third order result in a can be ob- 
served only for temperatures well below Ec/Wks even 
at ncx — Finally, Fig. |5| shows that the resumma- 
tion of the leading logarithmic terms (||) does not suffice 
to describe the behavior near riox — 5- Subleading loga- 
rithms are important to obtain quantitatively meaningful 
results. 
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FIG. 5. The average electron number (n) as a function 
of the dimensionless voltage n^^ is shown in first, second, 
and third order perturbation theory in a, and compared with 
QMC data. The result (8) is also shown as a dotted line. 

In summary, we have studied the breakdown of 
Coulomb charging effects with increasing tunnel conduc- 
tance. Precise Monte Carlo data for the effective charging 
energy as well as for the smeared staircase function were 
presented. Comparing analytical results with these data, 
the range of validity of expansions in powers of the tun- 
neling conductance was determined. It was found that 
presently attainable experimental results are covered by 
expansions up to third order in a. With increasing a, 
the effective charging energy i?* decreases rapidly. Since 
finite temperature corrections for large a are controlled 
by the dimensionless parameter /3£^* experiments 
with dimensionless conductance above a — 10, where 
perturbation theory begins to fail, require extremely low 
temperatures to see nonperturbative effects. We finally 
mention that the methods employed here can be extended 
to other systems displaying charging effects, such as the 
single electron transistor. 
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